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Summary of Observations 
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TABLE 1— Continued 
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Clean map. Array: BFHKtMNOPS 
0945+408 at 22.225 GHz 1995 Jan 04 
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Clean map. Array: BFHKtMNOPS 
1954+513 at 22.225 GHz 1995 Jan 04 
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Clean map. Array: BFHKLMNOPS 
2037 + 511 at 22.225 GHz 1995 Jan 04 



Clean map. Array: BFHKLMNPS 
2201+315 at 22.225 GHz 1995 Jan 04 
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VLBI Antennas 
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1993.18, Session 7:1995.01, Session 8:1995.34, Session 9:1996.34 
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TABLE 3 
VLBI Map Parameters 
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TABLE 4 
VLBI Gaussian Model Parameters 
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O A 

3.9 


A TO 

0.73 


O T O 

37.3 










A O "7 1 A A 1 

0.27 ± 0.01 


T AO 

7.03 


O A A O 

249.8 


A AT 

4.07 


A O O 

0.38 


1 E 1 

-15.1 








22 


r\ r> -i 1 A O 1 

0.61 ± 0.21 


0.0 


0.0 


0.35 


0.53 


12 


1.27 


1.15 






A A T 1 A 1 f 

0.47 ± 0.16 


2.5 


-109 


0.88 


1.0 


69 










A O A 1 All 

0.34 ± 0.11 


5.5 


-108 


1.4 


1.0 


88 






A r r o i o A o 

0552+398 


8.4 


/i at 1 A O O 

6.95 ± 0.23 


0.0 


0.0 


0.60 


0.82 


-7.3 


3.17 


1.59 




22 


E O O 1 A O A 

5.88 ± 0.34 


0.0 


0.0 


0.30 


0.77 


-6.0 


2.05 


1.12 


A/" 1 /I O i 1 A l~\ 

0642+449 


22 


O F* E 1 A O 7 

2.55 ± 0.27 


0.0 


0.0 


0.21 


0.71 


-91.6 


1.06 


1.27 


0736+017 


n o 

2z 


1 E O I A A T 

1.58 ± 0.47 


A A 

0.0 


A A 

0.0 


A AT 

0.07 


AO/ 1 

0.36 


OA 

29 


A OT 

0.87 


1 OT 

1.27 






A A E 1 A 1 /I 

0.45 ± 0.14 


0.52 


-78 


0.41 


0.48 


-19 










A 1 £? I A A E 

0.16 ± 0.05 


2.30 


-83 


3.15 


0.001 


48 






AO A /I i A AA 

0804+499 


O -1 

8.4 


A AT 1 A A/ 1 

0.97 ± 0.06 


A A 

0.0 


A A 

0.0 


A A A 

0.40 


A f A 

0.60 


T A A 

7.04 


1 CZ A 

1.60 


1.11 






A AA 1 A AE 

0.09 ± 0.05 


1 OO 

1.82 


TO A 

72.0 


A AO 

4.U2 


A 1 O 

0.13 


1 1 A 

11.0 








O O 

22 


1 *2 T I A O E 

1.67 ± 0.25 


A A 

0.0 


A A 

0.0 


A OT 

0.27 


A 1 A 

0.19 


1 o o 

16.8 


1 A A 

1.00 


1 A C\ 

1.40 


AOAA i e£?a 

0820+560 


22 


C\ A A 1 A AT 

0.44 ± 0.07 


A A 

0.0 


A A A 

0.00 


A O E 

0.25 


A OO 

0.23 


-67.5 










A 1 I A A O 

0.21 ± 0.03 


0.78 


105.6 


2.97 


0.13 


-42.1 






0836+710 


8.4 


1 A A 1 A A 1 

1.44 ± 0.01 


0.0 


0.0 


0.35 


0.46 


61.6 


2.76 


2.28 






A 17 1 A A 1 

0.47 ± 0.01 


2.33 


210.0 


1.65 


0.61 


14.5 








22 


1 TT I A 1 O 

1.77 ± 0.13 


0.0 


0.0 


0.18 


0.26 


32.1 


1.64 


1.07 


0917+449 


22 


1.31 ± 0.07 


0.00 


0.0 


0.20 


0.50 


7.4 


0.27 


0.21 






0.23 + 0.01 


1.13 


174.4 


0.99 


0.48 


14.6 






0945+408 


8.4 


1.75 + 0.25 


0.0 


0.0 


0.27 


0.38 


-74.6 


2.43 


1.43 




22 


1.57+0.03 


0.0 


0.0 


0.34 


0.59 


-20.5 


1.05 


1.07 




22 


1.09 + 0.05 


0.13 


105.0 


0.64 


0.25 


-82.8 


0.358 


0.304 


0955+476 


22 


1.01 + 0.05 


0.00 


0.0 


0.12 


0.89 


79.17 


0.742 


0.266 


1150+497 


22 


1.72 + 0.09 


0.00 


0.0 


0.19 


0.22 


19.9 


1.132 


0.101 






0.18 + 0.01 


0.61 


-155.5 


1.32 


0.50 


31.2 






1611+343 


22 


1.56 + 0.34 


0.0 


0.0 


0.24 


0.38 


1.4 


4.4 


1.33 






0.26 + 0.06 


0.45 


1.68 


0.33 


0.001 


5.3 







1 



TABLE 4— Continued 



Source 


Frequency 
(Ijrlz ) 


Flux Density 

( J y) 


r 

(mas) 


n 

"pos 

(degrees) 


a 

(mas) 


0/ a 


<P 

(degrees) 


A. r . 

(amplitude) 


A. r . 

(closure phase) 




22 


1.86 + 0.09 


0.08 


106.7 


0.22 


0.37 


-5.7 


0.309 


0.186 






0.57+0.03 


2.81 


171.1 


1.61 


0.49 


-38.7 










0.57+0.03 


0.15 


37.1 


6.86 


0.17 


50.0 










0.09 + 0.01 


0.86 


93.8 


1.19 


0.10 


-21.8 






1633+382 


8.4 


1.60+0.04 


0.0 


0.0 


0.52 


0.0 


84.6 


1.29 


1.20 






0.41 + 0.01 


0.84 


-75.7 


0.95 


0.89 


105.7 








22 


0.79 + 0.11 


0.0 


0.0 


0.24 


0.62 


37 


1.52 


1.29 






1.18 + 0.18 


2.6 


-32 


3.16 


0.40 


-32 






3C 345 


8.4 


9.41 + 0.09 


0.0 


0.0 


0.50 


0.12 


87.5 


3.56 


3.65 




22 


9.26 + 0.28 


0.0 


0.0 


0.33 


0.78 


2.9 


7.83 


1.35 






2.39 + 0.07 


0.29 


-87 


0.17 


0.59 


-6.4 






1954+513 


22 


0.72 + 0.11 


0.00 


0.0 


0.103 


0.10 


66.3 


0.438 


0.190 






0.41 + 0.06 


0.48 


-53.5 


0.36 


0.57 


-21.4 






2005+403 


22 


1.59 + 0.08 


0.0 


0.0 


0.80 


0.14 


83.0 


0.918 


0.498 






0.35 + 0.02 


0.51 


75.4 


0.83 


0.20 


60.2 






2037+511 


22 


1.98 + 0.10 


0.00 


0.00 


0.18 


0.52 


-6.5 


0.55 


0.331 






0.68 + 0.03 


0.75 


-157.1 


0.88 


0.37 


83.8 






2136+141 


8.4 


1.70 + 0.03 


0.0 


0.0 


0.14 


0.65 


-60.5 


4.63 


1.23 






0.45 + 0.01 


0.52 


-133 


1.96 


0.94 


-24.4 








22 


1.89 + 0.21 


0.0 


0.0 


0.33 


0.5 


-23 


1.4 


1.0 






0.49 + 0.05 


0.85 


-49 


5.7 


0.4 


-10 






2201+315 


22 


1.26 + 0.06 


0.0 


0.0 


0.23 


0.25 


31.3 


0.852 


0.457 






0.48 + 0.02 


1.87 


-142.9 


2.20 


0.28 


42.2 







2 




§ 10.5 



D) 
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TABLE 5 

Millimeter/ Submillimeter Observations with the JCMT c 



S\/"yi i T*r*iO 

OU Lll Lt 




V V dvclcllji 111 




Vpn-inhlp? d 

V cxl IclUlC 1 






w 


\ J y) 








(a J 


(A"\ 


(?) 


0736+017 


1993 Nov 21 


800 


1.13+0.08 


0.2 






1100 


1.27+0.08 


0.3 






1300 


1.48+0.11 








2000 


1.70+0.10 


0.3 


0917+449 


1993 Dec 08 


800 


0.84+0.10 


No 






1100 


0.95+0.10 


No 






1300 


1.01+0.10 








2000 


1.19+0.15 


No 


0923+392 


1993 Dec 08 


800 


0.48+0.06 


No 


(4C39.25) 




1100 


2.46+0.17 


No 






1300 


2.79+0.15 








2000 


3.62+0.45 


No 


0945+408 


1993 Dec 08 


800 


<0.39 


No 






1100 




No 






1300 


0.39 ±0.10 


No 






2000 


0.59 ±0.13 


No 


0954+556 


1993 Dec 08 


800 


<0.23 


No 






1100 


0.16±0.04 


No 






1300 


0.20±0.04 


No 






2000 


0.47±0.08 


No 


0955+476 


1993 Dec 08 


800 


0.21±0.07 


No 






1100 


0.23±0.05 


No 






1300 


0.40±0.05 


No 






2000 


0.44±0.08 


No 


1641+399 


1993 Sep 14 


800 


1.35± 0.12 


0.4 


(3C 345) 




1100 


1.85± 0.07 


0.3 






1300 


2.24± 0.10 








2000 


2.91±0.15 


0.3 


2005+403 


1993 Sep 14 


800 


0.23±0.07 








1100 


0.29±0.05 


0.3 






1300 


0.47±0.08 


No 






2000 


0.70 ±0.14 


No 


2201+315 


1993 Sep 14 


1100 


0.40 ±0.07 


0.3 






1300 


0.41 ±0.06 


0.3 






2000 


0.60 ±0.13 


0.3 



c see Bloom et al. 1994 for data from earlier epochs 

d F ma x-F mi „ ."jsjq"^ if variability is not significant, as defined in text 



1 



log F„ X-RAY FLUX DENSITY (/zJy) 




TABLE 6 
Infrared Observations with UKIRT 



Source Date Wavelength Flux Density aiR Instrument 

(//m) (mJy) 







(3) 










y°) 


0804+499 


1993 Nov 19 


1 


o r 


1 


70+0 


00 


-0.01 


UKT9 






1 


c r 


1 


81+0 


07 










2 


20 


1 


69+0 


04 








1994 Mar 5 


1 


25 





33+0 


01 


0.06 


UKT9 






1 


05 





41+0 


03 










2 


20 





37+0 


03 






0820+560 


1993 Dec 09 


1 


25 





30+0 


02 


0.63 


UKT9 






1 


05 





42+0 


04 










2 


20 





43+0 


03 






0923+392 


1993 Nov 19 


1 


25 


1 


59+0 


03 


0.38 


UKT9 


(4C 39.25) 




1 


05 


1 


43+0 


03 










2 


20 


1 


97+0 


07 






0945+408 


1994 May 13 


1 


o r 

25 





33+0 


00 


0.90 


IRCAM 






1 


05 





45+0 


09 










2 


20 





55+0 


11 








1994 May 17 


1 


25 





30+0 


00 


0.97 


IRCAM 






1 


05 





39+0 


08 










2 


20 





52+0 


10 






0954+556 


1993 Dec 09 


1 


25 





77+0 


02 


0.94 


UKT9 






1 


05 





96+0 


03 










2 


20 


1 


31+0 


04 






0955+476 


1993 Dec 09 


1 


25 





71+0 


03 


1.61 


UKT9 






1 


05 


1 


05+0 


05 










2 


20 


1 


76+0 


07 






1641+399 


1993 Sep 15 


1 


25 


1 


03+0 


01 


1.51 


UKT9 


(3C 345) 




1 


05 


1 


48+0 


06 










2 


20 


2 


42+0 


06 






2005+403 


1993 Sep 15 


1 


25 





29+0 


02 


0.66 


UKT9 






1 


05 





33+0 


02 










2 


20 





42+0 


04 






2201+315 


1993 Sep 15 


1 


25 


3 


99+0 


03 


1.52 


UKT9 






1 


05 


5 


48+0 


03 










2 


20 


9 


39+0 


05 







a see Bloom et al. 1994 for data from earlier epochs 



log L x (erg/sec/Hz) 
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TABLE 7 
X-ray Spectral Parameters 



Source 




Date 


Counts 
(0.1-2.5 keV) 


Obs Time 
(sec) 


log N H ,gal 

(atoms cm -2 ) 


log N H ,bcst 

(atoms cm~ 2 ) 




a 


FlkeV 
(A<Jy) 


x 2 


0016- 


-731 




Jul 29 


40+7 




21 


38 




\ 


n 

u 


n OR + o m 

U.UU —1— U.U1 


fi Q 

U . L? 


0133- 


h476 


1 QQ9 


A ncr 19 
r\U^ 1Z 


oat; _i_ i 7 
zoo m i ( 




21 


37 


on Q9+008 
zu.»z_q 12 


n 


fifi +0.64 


n 97+0.10 

u - z '-0.04 


1 _L 

1-i 


0212- 


h735 a 


1 QQ1 




oris -1-18 


791 


21 


39 








n 9n+ 02 

U./U_ .03 


Q 


NRAO140 b 


1992 


Aug 8 


921+31 


4039 


21 


41 


21.46 ± 0.02 


1 


01-0.17 


- r „-l_l.ll 

i-SSIo.ei 


33 


3C lll b 


1991 


Mar 8 


265 + 17 


838 


22 


08 




2 


or+0.65 
AO -0A0 


12.3j£» 


7 






1993 


Feb 13 


840 ± 29 


2322 


22 


08 


21.98 + 0.02 


1 


48 + 0.20 


8.4™ 


16 


3C 120 


1991 


Feb 20 


3642 ± 62 


1724 


21 


03 


21.54 + 0.02 


2 


-i o+0.09 
io -0.11 


18.4i^ 


36 


0552- 


-398 


1993 


Mar 13 


185+14 


4294 


21 


28 







nq +0.36 
uy -0.16 


0.21 + 0.03 


12 


0642- 


-449 


1991 


Mar 13 


73+10 


3108 


20 


79 







7O+0.86 
' y -0.59 


0.08t°;£ 


20 


0736- 


-017 


1991 


Apr 3 


234+16 


1936 


21 


08 




1 


0/1 +0.36 
z ^-0.14 


U.OU_ 05 


21 






1992 


Apr 6 


396+20 


3296 


21 


08 




1 


4C+0.25 
^ d -0.45 


n 47+0 06 


24 


0804- 


h499 


1992 


Oct 27 


144 ± 12 


2999 


21 


32 




1 


nn+0.71 
uy -0.79 


24+ - 11 

u - z ^-0.07 


44 


0836- 


-710 a 


1992 


Mar 3 


5365 ± 75 


6990 


20 


46 


20.54+ 0.05 





57 + 0.11 


2.37+0.10 


41 






1992 


Nov 2 


1904 ± 44 


5030 


20 


46 


20.56+ 0.10 





t- 7 +0.19 
°' -0.18 


-, 91 +0.08 
i - zl -0.09 


15 


0917+449 


1993 


Oct 27-Nov 3 


639 ± 28 


3225 


20 


55 


20.53+ 0.10 





y *-0.24 


n 4Q+0- 25 


40 



1 



TABLE 7— Continued 



Source 


Date 


Counts 

\\J.L Z.O KC V J 


Obs Time 

(sec) 


log N H ,gal 

f of ATTIC /~*VVI ^\ 

(^aioms cm j 


log N H ,bcst 

^aioms cm j 


a 




FlkeV 


x 2 


D. 


nnoQ i OOOC 


iyyi Apr io 


1 41 




on oo 


on on _i_ n on 
zU.zU ± U.zU 


a i 7+0.38 
L17 -0.34 


n 
U 


ko j_ n nf; 
.Dz ± U.Ub 








ryyo Apr zz 


1 Q CCf? 

loOD 




on oo. 
zU.zo 


on qn+0.16 


1 qo 4- n Qn 
l.oz ± U.oU 


n 
u 


gq+0.07 
by -0.06 






0y45+408 


i nno at^,. 1 n 

iyy3 INOV 1U 


1 on i 1 
189 ± lb 


AC77 

25 (V 


on i £; 
zU.lb 




1 n«+0.22 
1.06I 024 


n 



nn _i_ n no 
.09 ± O.Uz 


a a 

4b 


O 1 

31 


0y54+55b 


1 nno at^,. o 

iyy3 INOV 8 


ono _L O A 

3y3 ± 24 


a ono 
4293 


1 n n c 

iy.95 


19.92_ 040 


-i o.n+0- 45 

1 - ,3u -0.30 


n 



no _i_ n no 
.08 ± 0.02 


c o 

53 


O 1 

31 


U955+4/D 


iyy3 INOV 13 


o on _i_ on 
33y ± 20 


O A C A 

345U 


on m 
20.01 




n i n+0.15 
U ' ly -0.26 


n 



91 +0.05 
• zl -0.04 


A K 

45 


O 1 

31 


1 1 rA i o 1 o a 

1150+812 


i nno A , n 

iyyi Apr r 
1992 May 7 


c o _i_ n 

58 ± y 

172+14 


onon 
293U 

5650 


on a c 
20.65 

20.65 




n crc+0.21 
U.OO_ 25 

(Combined) 


n 



nn+0.09 

uy -o.io 


1 o 

13 


5 




1992 Oct 29 


241 + 17 


9500 


20.65 




41 + - 19 





QO+0.10 
uy -0.09 


5.4 


4 


1611+343 


1992 Sep 30 


418 + 21 


3651 


20.45 


20.32^°,;°! 


1.34+0.30 





1 r+0.03 

' lo -0.01 


80 


30 




1993 Jan 17 


332+ 18 


3153 


20.45 


20.30±°,; 25 ' 


0.65 + 0.35 





94+0.04 
• z ^-0.03 


lb 


24 


1633+382 


1992 Aug 20 


1208 + 38 


5781 


20.58 


20.18±8;l| 


76 +0 17 

u - '°-0.20 





4 i+0.04 
' 41 -0.03 


41 


30 


1638+398 


1992 Sep 3 


64+9 


2857 


20.49 




1.89 + 0.30 





.04+0.01 


15 


15 


3C 345 


1990 Jun 19 


1193 + 43 


4086 


20.56 


iq 71 +o.20 


n 77+0.13 

U.I <_ .17 





qn+0.04 
oy -0.03 


54 


30 




1992 Jun 16 


1014 + 40 


3899 


20.56 


20 n+ 014 


u - yi -0.25 





.51 + 0.05 


41 


30 


2136+141 


1991 May 9 


147+16 


4076 


20.93 




O0+ - 50 





.12 + 0.02 


22 


20 


2201+315 


1993 May 12 


553 ± 25 


3329 


21.09 


21 29 +nn 


2 00+ - 90 
z.uu_ 70 


3 


79+0.62 
''°-0.58 


14 


32 



a from Brunncr et al. 

b Nn is atomic and molecular 

c from Zhang et al. 



2 



TABLE 8 
Spectral Dissection Parameters 



Source 


v n 




1 m 




Otthick 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0016+731 


60 


54 


1.94 


1.36 


0.33 


0133+476 


42 


37 


3.80 


0.65 


1.13 


0212+735 


64 


58 


1.99 


1.57 


0.0 


NRAO 140 


78 


69 


2.92 


1.35 


0.38 


3C 111 


79 


73 


4.28 


0.75 


0.81 


3C 120 


67 


67 


2.63 


0.80 


0.41 


0552+398 


39 


47 


4.35 


1.06 


0.0 


0642+449 


75 


65 


4.09 


1.81 


0.0 


0736+017 


58 


78 


1.95 


0.84 


0.12 


0804+499 


45 


57 


2.15 


0.63 


0.56 


0836+710 


27 


33 


1.34 


0.90 


0.21 


0923+392 




8 


4.5 


0.52 




0945+408 


33 


33 


1.73 


0.72 


0.0 


0954+556 


19 


28 


0.88 


0.65 


0.28 


1611+343 


47 


51 


3.18 


0.65 


0.88 


1633+382 


74 


94 


1.82 


0.60 


0.56 


3C 345 


48 


53 


8.98 


0.60 


0.90 


1954+513 


27 


38 


2.52 


0.65 


0.35 


2005+403 


37 


26 


4.18 


0.92 


1.58 


2037+511 


31 


50 


3.88 


0.63 


0.11 


2136+141 


19 


15 


1.63 


0.77 


0.27 


2201+315 


26 


21 


4.32 


0.70 


1.47 



1 



TABLE 9 
Compton Calculation Parameters 



Source 


z 


% 


± m 




9 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


0016+731 


1.76 


54 


1.94 


1.36 


0.67 


5.1 x 10~ 8 


0.06 


0.13 


0133+476 


0.86 


37 


3.80 


0.65 


£0.os 


£ 2.6 x 10 5 


0.27 


14 


0212+735 


2.37 


58 


1.99 


1.57 


0.36 


7.3 x 10~ 7 
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ABSTRACT 

We present multiwaveband observations of a well-selected sample of 28 
quasars and two radio galaxies with flat radio spectra and strong millimeter-wave 
emission (referred to here as FSRQ's). The data are analyzed to determine the 
radio to infrared and X-ray to 7-ray properties of FSRQ's and the relationships 
between them. Specifically, the synchrotron self-Compton (SSC) process is 
examined as a likely common radiation mechanism. For most sources, the broad 
band spectra are still incomplete, especially in the far-infrared and ultraviolet 
range. Therefore precise analysis, such as model-fitting of spectra is not usually 
possible. To compensate partially for this, we have taken a statistical approach, 
and examine the relationship between high and low energy emission by using 
the data set for the entire sample. 

We use very long baseline interferometry (VLBI) at 8.4 and 22 GHz - 
higher frequencies than those of previous surveys — in conjunction with 
nearly simultaneous radio to submillimeter-wave observations to determine the 
parameters of the synchrotron spectrum and to examine the compact angular 
structure of a subset of sources from our sample. These parameters are used 
to predict the SSC X-ray flux densities. Seven of thirty sources have predicted 
self-Compton X-ray flux densities well above the observed flux densities 
obtained with the ROSAT satellite unless one assumes that the radiating plasma 
experiences bulk relativistic motion directed toward the observer's line of sight. 
Three of these seven sources are detected at 7-ray frequencies. Model spectra 
show that the X-rays are consistent with the first order SSC process, with the 
simultaneous multiwaveband spectrum of the quasar 0836+710 obtained in 
1992 March being very well fit by SSC emission from a uniform, relativistically 
moving source. The 7-rays are not produced via second order self-Compton 
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scattering, but rather by either first order self-Compton scattering or some other 
process. 

A comparison of the ROSAT X-ray flux densities and those obtained earlier 
with the Einstein Observatory show that several FSRQ's are X-ray variables on 
timescales of about a decade. Several sources that were observed more than 
once with ROSAT also show variability on timescales of 1-2 yr, with the X-ray 
variability in these cases often associated with millimeter-wave variability and 
lower VLBI core-to-jet flux ratios. Detections at 7-ray energies also appear to 
be related to increases in the radio to millimeter-wave flux densities. 

Statistical analysis shows that the millimeter-wave and X-ray luminosities 
for the sample are strongly correlated, with a linear regression slope ~ 0.6. The 
peak in the distribution of X-ray to millimeter spectral indices also indicates 
a strong connection between millimeter-wave and X-ray emission. Particularly 
interesting is a correlation between X-ray to millimeter spectral index and 
fraction of flux density contained in the VLBI core. This tendency toward 
higher X-ray fluxes from sources with stronger jet emission implies that the 
knots in the jet are prominent sources of X-rays. 
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1. Introduction 

It is clear that in quasars there is a strong relationship between the nonthermal emission 
at low frequencies and that at high frequencies. Quasars are generally strong X-ray emitters 
(e.g., |Tananbaum et al. 1979| ; |Ku et al. 1980| ), and a number have been found to be strong, 



hard (> 100 MeV) 7-ray sources QFichtel et al. 1994| [Thompson et al. 1995| ; [Thompson et 



al. 1996|) . Quasars with flat radio spectra — which indicate compact, milliarcsecond-scale 



structure — and strong millimeter-wave emission (here abbreviated as FSRQ's) have been 
demonstrated to be more luminous in soft X-rays (for a given optical luminosity) than 
either their radio-quiet counterparts or radio sources with the same radio luminosity but 
steeper radio spectral index ( Pwen et al. 198"!] ; |Zamorani et al. 198 1| ; |Ledden et al. 1985 ; 



[Kembhavi et al. 1986| [Worrall et al. 19"87|) . |Browne et al. 1987] have also shown that for a 



given radio luminosity in extended (arcseconds or greater) structure, the X-ray luminosity 
tends to be more than an order of magnitude larger for sources with higher compact radio 
luminosity. Furthermore, except for the Large Magellanic Cloud, bright extragalactic hard 
7-ray sources are all FSRQ's (or BL Lac objects, which have similar radio-infrared spectra). 

These strong correlations between radio and high-energy emission are readily explained 
if the nonthermal emission at all wavebands is produced via synchrotron self-Compton 
(SSC) scattering, as has been suggested by many authors (e.g., [Jones et al. 1974a|). This 



is, however, not the only possibility, since other radiative processes are likely to operate in 
an environment that produces the highly relativistic electrons and strong magnetic field 
required for SSC emission (eg., permer et al. 1993| ). 



If the SSC process dominates the nonthermal emission, one would expect a correlation 
between X-ray and 7-ray fluxes and VLBI "core" fluxes. However, a statistical study of 
the (non-simultaneous) VLBI core emission and X-ray emission from compact radio loud 
quasars and active galaxies ( [Bloom fc Marscher 1991] ) did not conclusively demonstrate a 
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strong correlation. Though Zhou et al. (1997) do find a strong relationship between high 
energy 7-ray fluxes and VLBI fluxes. Nevertheless, detailed studies of individual sources, 
both with VLBI and X-ray data, show that the SSC process is a plausible explanation for 
the X-ray emission ( |Marscher 1988| ; [Unwin et al. 1994| ; |Unwin et al. 1997| |Eckart et al.| 



1986| ; |Eckart et al. 1987|) . The results of these investigations have motivated us to conduct a 



more detailed study of FSRQ's and the relationship between the radio-infrared and the high 
energy radiation. The present study not only explores the statistical relationships between 
high and low energy emission, but also compares theoretical SSC spectra to multiwaveband 
data from a set of contemporaneous measurements. Observations over a large range in 
frequency allow for more accurate determinations of spectral index, spectral peak frequency 
and flux density, and consequently a more accurate prediction of the Compton flux densities. 

This study, which includes high-frequency (8.4 and 22 GHz) VLBI imaging, as well 
as radio through 7-ray (not all wavebands for all sources) flux density measurements, 
constitutes the most comprehensive set to date of multiwaveband observations of bright 
FSRQ's. Below we summarize the observations and data analysis. We present our own 
VLBI, X-ray, submillimeter-wave (submm), and near- infrared (IR) data. A more extensive 
presentation of previous radio-IR observations used in this work is given by [Bloom et 



|al. 1994 We follow with explanations of the procedure for spectral deconvolution and 



determination of spectral peak frequency and peak flux density, SSC calculations for specific 
sources, and comparisons of model spectra against the data. In the final sections we present 
a statistical analysis of the sample and a general discussion of the results. 



2. Observations and Data Analysis 
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2.1. Source Selection 

The sample is taken from an essentially complete list of radio sources stronger than 1 
Jy at 5 GHz, with spectral indices between 1.5 and 4.9 GHz flatter than a = 0.5 (F v oc v~ a ) 
( [Steppe et al. 19881) . Quasars and galaxies with 90 GHz flux densities exceeding 0.5 Jy, 
declinations > 0°, and measured redshifts exceeding 0.03 were selected from Table III of 
(Bteppe et al. 1988| to generate a preliminary sample. This resulted in a sample that had a 



preponderance of sources with redshifts between 0.4 and 2.25. In order to pare the sample 
down to a more reasonable size, within the redshift range of 0.4 to 2.25 only sources with 
declinations > 30° were included in the final sample. BL Lac objects, which form a class 
distinct from quasars with respect to their soft X-ray emission ( [Worrall fc Wilkes 1990|) are 



excluded, as are optically unidentified objects. The final sample contains 30 sources — 28 
quasars and two radio galaxies — and is well distributed in redshift. (One source that would 
have been in the final sample, 3C 216, was dropped from the study due to its proximity 
to the Sun during the time of the submillimeter observations.) Because of telescope time 
allocation limitations, it was not possible to observe the entire sample at all wavelengths. 
Approximately one half of this sample was observed from radio to X-ray wavelengths. This 
subsample was essentially determined by which sources were observed (in priority order) 
with ROSAT prior to the shutdown of the PSPC detector. This might then have caused 
a selection of objects that were a priori deemed to be more interesting than others. Both 
radio galaxies, 3C 111 and 3C 120, are included in the subsample. The list of sources, along 
with a summary of the observations at all wavebands, is presented in Table 1. 



2.2. VLBI Data 

A subsample of sources has been observed with VLBI at frequencies of 8.4 and 22 
GHz. For this work, seven observing sessions were conducted with either the Global VLBI 



- 7- 



Network or the Very Long Baseline Array Q( VLB A; after 1994). The stations of each array 
(defined by the epoch of the observations) are listed in Table 2. Sessions 1-6 used the Mark 
II VLBI recording system with a 1.8 MHz bandwidth, while sessions 7-9 used the VLBA 
recording system with 16 and 32 MHz bandwidth respectively. Due to the large sample 
size, we observed in "snapshot" mode with each source typically observed for 4 to 5 scans of 
one-half hour each during the allotted observing time. Typically three to six sources were 
observed during each VLBI session. These scan times were chosen to optimize u-v coverage. 
The video tapes from the individual antennas were correlated using the JPL/Caltech Block 
II correlator for sessions 1-6 and the VLBA correlator for sessions 7-9. All post-processing 
was performed at Boston University using the NRAO AIPS, the CalTech VLB software 
package, and the difference mapping software Difmap. The phase time and frequency 
derivatives (fringe delay and fringe rate) were determined using the global fringe fitting 
algorithm FRING in AIPS ( |3chwab fc Cotton 1983| ). The amplitudes were calibrated in 
the usual manner as described by Cohen et al (1975). Inconsistencies in the calibration 
of flux densities were corrected by using u-v crossing points and then finding the best-fit 
flux density using a x 2 test. A hybrid mapping procedure was used to obtain images of the 
sources (see, e.g., [Pearson fc Readhead 198-31 ). Hogbom's CLEAN algorithm was used to 
deconvolve the dirty beam from the dirty image to determine the true source brightness 
distribution. A Gaussian restoring beam with dimensions equal to those of the central 
portion of the dirty beam was convolved with the CLEAN components to produce a final 
map. 

Table 3 summarizes the parameters of the final image for each source at each frequency: 

x The VLBA is an instrument of the National Radio Astronomy Observatory, a facility 
of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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elliptical Gaussian restoring beam (abeam, Pbeam), the position angle of the beam (6*beam, as 
measured north through east), the peak flux density per unit beam area, and the contour 
levels. Table 4 summarizes the Gaussian model parameters used in conjunction with the 
hybrid mapping procedures, as described above. These are from the model fits to the 
self-calibrated data. As with the Pearson & Readhead (1988) survey, the limited data only 
weakly constrain some model components. In Tables 3 and 4, r is the angular distance from 
the map center, ^pos is the angular position of the component relative to the (arbitrarily 
placed) origin (measured from north through east), a and b are the FWHM major and 
minor axes of the assumed elliptical Gaussian brightness distribution of the component, 
and <p is the position angle of the major axis. The uncertainties in the flux densities given 
are the estimated calibration errors (as opposed to the uncertainty in the model, which 
is usually larger). Though we have not, in general, estimated the uncertainties for every 
parameter of each component, we have performed an uncertainty analysis on the angular 
sizes (characterized as major and minor axes of an elliptical Gaussian here) of some sources. 
This will be helpful later (§3.2) in determining the uncertainty of the predicted Compton 
flux, which is highly sensitive to the angular size value. In determining the magnitude of 
the 1-a uncertainities, we use the technique of Biretta et al ( 1986) Due to the large range in 
the data quality of our snaphot images, there is also a fairly large range in the uncertainties 
of derived angular sizes. For the 4 objects observed with the VLBA and also later used for 
the Compton calculations, the uncertainties are 13/for 0945+408, 0955+476, 1611+343, 
and 2201+315 respectively; however, for sources observed earlier with the Global VLBI 
Network, the typical model uncertainties are several times larger. Figure 1 displays the 
VLBI images at 8.4 and 22.2 GHz, in order of increasing right ascension for sources mapped 
with Global VLBI. Figure 2 displays the maps from the VLBA. Note the differences in 
angular scale, especially between most of the 22 GHz and 8.4 GHz maps. 
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2.3. Submillimeter and Infrared Data 

A detailed description of the data analysis at submillimeter and infrared wavebands is 
given in Bloom et al (1994). In Tables 5 and 6 we present data more recently acquired with 
the JCMT and UKIRT. The last column in Table 5 is a measure of variability amplitude, 
as determined by the parameter ^ max ^gjg; ■ The two measured flux densities for this 
calculation, extracted from the data presented here or in Bloom et al (1994), are from the 
two observations that are closest in time. This measurement of the variability amplitude is 
considered to be significant if there is at least a 2-cr separation between the error bars of the 
highest and lowest point. We have also examined the possibility of variability in the older 
data (between 1991 and early 1993). The sources found to be variable are discussed in §2.6. 



2.4. X-ray Data 

The Position Sensitive Proportional Counter (PSPC) was used in conjunction with 
the Woltjer Type I nested mirror system onboard the Roentgensattelit (ROSAT) to obtain 
fluxes over the energy range 0.1-2.5 keV. After initial processing at Goddard Space Flight 
Center (GSFC), the X-ray data were extracted from tape using the Post Reduction Off-line 
Software (PROS) routines operating under the NOAO IRAF data analysis package. Before 
conducting a spectral analysis, the background counts were determined for each source by 
calculating the number of photons within concentric annuli, centered on the source. The 
precise sizes of the annuli (typically 2' inner and 3' outer) were determined for each source 
by displaying the photons/pixel with SAOIMAGE (see [Wilkes 1992 ). All of the objects in 



this sample are unresolved with the PSPC, which has a spatial resolution of about 15". 
Therefore, spatial analysis of the photon counts (i.e., mapping the source) is not possible. 

The spectral analysis was conducted using the standard routines within PROS. This 
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consisted of modeling the data in terms of a power-law spectrum over the energy range 
0.1-2.5 keV, with absorption at low energies due to intervening gas, parameterized as the 
column density of hydrogen atoms under the assumption of cosmic abundances. The X-ray 
photon spectrum is represented by the following equation: 

= KE -{a+i) e -N H a(E) photons s^cnrteV- 1 . (2.4.1) 

dE 

The predicted spectrum (photon counts per energy bin) using this model was then compared 
to the observed spectrum using a \ 2 test to determine "goodness-of-fit." The model 
parameters that minimize \ 2 were adopted to calculate the total flux over a specified energy 
band. The default values for the spectral response matrix and effective area of the detector 
for routines in PROS 2.0 were used to determine the source spectrum from the raw photon 
spectrum. Some care has to be taken, as x 2 can be a minimum for parameter values that are 
physically unrealistic, such as an extremely high column density, or a photon index that is 
much too steep. Uncertainties in model parameters and derived flux values were determined 
by increasing the value of the parameter by the amount that increases \ 2 to the desired 
confidence level. Alternatively, in some cases, the normalization parameter, K, of equation 
(2.4.1) was plotted against spectral index, and the uncertainty in this parameter was taken 
to be the range of the 68% confidence contour. The flux density was then calculated by 
using equation (2.4.1) corrected for absorption, and then converting to Janskys. 

Table 7 summarizes the ROSAT X-ray observations. Column (5) is the logarithm of 
the Galactic hydrogen column density, iV# gal . Unless otherwise noted, this is the value 
from the Stark et al. (1992) 21-cm emission line observations. However, for two sources, 3C 
111 and NRAO 140, we have also added the column density (in atoms cm -2 ) of molecular 



hydrogen as inferred from observations of Galactic CO ( Bania et al 1991 ). Column (6) 
gives the best fit value of Nh as determined from the x 2 minimization routine. Column 
(7) is the best-fit spectral index; if no uncertainties are given, then it is an assumed value. 
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Column (8) is the 1 keV flux density. Columns (9) and (10) give the x 2 an d degrees of 
freedom (D.O.F., one less than the number of parameters varied subtracted from the total 
number of energy channels), respectively. When Nh is determined from a best fit, then 
the uncertainites quoted in the spectral index and flux densities correspond to Xmin +3.53 
(68% confidence for three interesting parameters). If N H is assumed to be the Galactic 
value, then the subsequent uncertainties correspond to Xmin + 2-3 (68% confidence for two 
interesting parameters). If both a and Nh are assumed then the uncertainties are for 
Xmin +1-0 (68% confidence for one interesting parameter). Generally, the fit is considered 
reasonable if x 2 ^ D.O.F. Note that by this criterion two fits are not particularly good (see 
Table 7). The fit for the first observation of 1150+812 is poor probably because the fit was 
performed on the combined counts from two scans that were a month apart, each having a 
low number of individual counts. Also, the actual Nh value could be smaller or larger than 
the assumed Galactic value (see below for a discussion of this point). It is also possible that 
a single power-law model is not adequate for describing this source. The first observation 
of 1611+343 is poorly fit. Probably only the last reason is applicable here, since there are 
enough counts to determine the spectrum rather well, and the best-fit N H is fairly close to 
the Galactic value. 

In many cases the best-fit values for N H are either considerably smaller or larger than 
the measured Galactic N H . There are several likely explanations for this. If the best-fit 
value is larger, as is the case for 3C 120 and 0736+017, then there could be significant 
photoelectric absorption within the source. There could also be significant amounts of as yet 
unobserved Galactic molecular gas in the direction of the FSRQ. Also possible is a spectral 
flattening at low X-ray energies. A closer look at the photon spectrum of 3C 120 shows 
that the peak is poorly fit by the model, indicating that a single power-law is probably 
not the best model for this source. It is also possible that the Galactic column density of 
hydrogen is higher (or lower) in the immediate vicinity of the FSRQ as compared to the 
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average value over the much larger regions measured with the 2° FWHM beam of the Stark 
et a/.(1992) survey. This seems to be the case for 3C 345, for which Nh measured with 
the smaller beam in the Elvis et al. (1989) study is 7.5 x 10 19 , which is certainly closer to 
the best-fit value of Nh for the 1990 observation than is the value Af# gal = 3.6 x 10 20 cm -2 
derived from Stark et al (1992). These possibilities have been discussed in more detail by 
Wilkes & Elvis (1987) and Worrall & Wilkes (1990), among others. However, the higher 
resolution HI observations of Elvis et al. (1989) show that A# gal is smaller than derived 
from Stark et al (1992) for 0736+017 and only 2% larger for 3C120. It is still possible 
that Nh s13 \ in the immediate direction of the quasar was larger at the time of the X-ray 
observations, since temporal variability of N HgSk \ cannot be excluded (see, e.g., [Fichtel et 



al. 1994 ). For 0736+017, the inferred column density from CO observations ( [Liszt et~aT 



1993| ) is 5.6 (±0.8) x 10 20 cm 2 , which can account for the difference. This also results 



in a much steeper derived value of the spectral index, which is more in agreement with 
other low-redshift active galaxies and quasars (see discussion below). If the best-fit value 
of Nh is smaller than the Galactic value, as is the case for 0133+476, 3C 111, 1633+382, 
and 3C 345, then it is possible that these sources are not adequately described by a single 
power-law. The most likely possibility is a soft X-ray excess, perhaps from thermal black 
body X-ray emission (from an accretion disk) with a steep slope. 

An analysis for some sources in this sample has appeared in other summaries of X-ray 
emission from blazars (see |Comastri et al. 1997[ |Sambruna 1997] ) . Most of the parameter 
determinations (such as spectral index and column density of hydrogen) are very similar. 
Our flux density values do occasionally differ, mainly because we base our flux values on 
best fit parameters, whereas Sambruna (1997) use fixed Af# gal to calculate flux densities. In 
addition, their source/background counts determination differs from ours slightly. The first 
observation of 1611+343 is the only case for which the flux densities differ by more than 5% 
(we present a value which is 50% lower). 
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2.5. Gamma- Ray Data 

Several sources in our sample were detected with the EGRET instrument on the 
Compton Gamma-Ray Observatory. The fluxes for EGRET detected objects can be found 
in Thompson et al. (1995) and upper limits for some sources in ths sample can be found 
in Fichtel et al. (1994). The dates listed in Table 1 are for the EGRET Viewing Periods, 
covering the area of the source, which were closest in time to the other mult i- wavelength 
observations. Calibration of the instrument and determination of fluxes are discussed in 
Hartman et al. (1992), Mattox et al. (1993), and Mattox et al. (1996), respectively. 



2.6. Summary of Source Properties 

Due to incomplete time coverage, it is not possible to analyze thoroughly the variability 
properties of each source. However, we can discuss some general results. Of the 17 sources 
with repeated millimeter/submillimeter wave observations reported here and in Bloom et 
al. (1994), ten (0133+476, NRAO 140, 3C 111, 0642+449, 0736+017, 0836+710, 3C 345, 
2005+403, 2037+511, and 2201+315) are significantly variable at one or more wavelengths. 
Six of these (0133+476, NRAO 140, 0736+017, 0836+710, 3C 345, and 2201+315) are 
significantly variable at all measured wavelengths. The timescales are roughly 7 to 14 
months; no source has a variability amplitude parameter that exceeds 0.5. 

With the X-ray observations presented earlier, combined with observations by previous 
spacecraft missions (HEAO A-l, Einstein, EXOSAT, and Ginga), we have determined 
the crude X-ray variability characteristics of the sources in the sample, using the same 
criterion for significance as for the millimeter observations. We note that in some cases it 
is difficult to establish variability, especially if the observations being compared were made 
with different instruments or very different spectral fit parameters (as mentioned above). 
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The sources appear to fall into several general categories ( Malaguti et al. 1994 ; Bloom 



1994| for references to the published X-ray data for each source). Of the 24 sources in 



our sample with ROSAT X-ray observations, seven (0133+476, 0552+398, 0820+560, 
0917+499, 0954+556, 0955+476, and 1638+398) have indeterminate X-ray variability 
due to insufficient observations. For an additional seven sources (0016+731, 0642+449, 
0736+017, 0804+499, 1150+812, 2136+141, and 2201+315) variability was not detected, 
though some of these are low amplitude variables at mm wavelengths (e.g., 0736+017; 
see [lbrnikoski et al. 1994 ). Two sources, 0945+408 and 1633+382, are observed to be 
significantly variable in X-rays (and, for 1633+382, 7-rays; see [Thompson et al. 1995| ) only, 
whereas three additional sources (NRAO 140, 0836+710 and 4C 39.25) show possible X-ray 
variability that is temporally correlated with radio or mm-wave variability. With each of 
these three sources, the mm flux varied by the same factor as the X-ray flux (see |Marscher 



19881 for NRAO 140; this work and Bloom et al. (1994) for 0836+710 and 4C 39.25). These 



last five sources mentioned were all variable with constant X-ray spectral index, which 
is typical of nonthermal processes. Five additional sources (0212+735, 3C 111, 3C 120, 
1611+343, and 3C 345) had complex variability properties, such as variation in spectral 
index as well as intensity. All of the sources with ^ 80 % of the VLBI flux density in the 
core at 8.4 GHz are X-ray variable, falling in one of the last three categories. The possible 
relationship between X-ray and milliarcsecond-scale structure is discussed in §4.2 . 



3. Spectral Analysis 

3.1. Spectral Deconvolution 

In this section, the radio to submillimeter and infrared spectra of Bloom et al. (1994) 
are used to determine the synchrotron spectral turnover frequencies for each source in the 
sample with sufficient data. In addition, for sources observed with both ROSAT and VLBI, 
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X-ray measurements are compared with the theoretical spectra of synchrotron self-Compton 
(SSC) models generated by the techniques described below. For sources with published 
EGRET 7-ray data, those data are also included for further comparison with the theoretical 
spectra. The total flux density spectra presented in Bloom et al. (1994) are dissected into 
the spectra of individual VLBI scale components using the total flux density data and flux 
densities from the best-fit models to the 8.4 and 22 GHz VLBI visibilities. The emphasis 
of the spectral decomposition is on finding the spectral turnover frequency and flux density 
of the bright "VLBI core" component usually located at one end of a source. Due to the 
limited data and the weakness of the secondary components we do not attempt to fit 
the spectra of any non-core "knots" (except for the peculiar source 4C 39.25, which was 
analyzed by [Zhang fc Marscher 1994| as part of a separate study) . 

In Figure 3 the spectral dissection for a representative source is shown. The turnover 
frequency v m is determined from the data by first calculating the frequency v n at which the 
optically thin and optically thick spectra (of slope ct mm and ctthick, respectively) intersect 
( Pones et al. 1974b| ; [Marscher 1977|) . In some cases this is not straightforward, since 



the VLBI spectra are complex (e.g., 3C 120 and 1633+382). For such sources, there is 
significant blending of components at 8.4 GHz that are separated in the higher resolution 22 
GHz maps. There are also additional components not visible in the 22 GHz images because 
of over-resolution and/or steep spectra. There are insufficient data to separate these 
components properly, thus only the 22 GHz core flux is used with the total flux density 
data to determine the spectral turnover. In sources with no VLBI data, best estimates are 
made using the total flux density data. We also note that in some cases, the VLA and/or 
VLBI data points exceed the values for the single-dish measurements (Figure 3). Since it is 
obviously not possible for the small scale components to have greater flux density than the 
entire source seen by a single dish, these differences in flux are attributable to differences 
in calibration. Since all of the sources have falling millimeter and submillimeter spectra, 
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there is a firm upper limit of 100-150 GHz to the turnover frequency. All sources in the 
sample have optically thick spectral indices below the turnover that are much flatter than 
the value of 2.5 expected from a homogeneous source. These flatter indices are indicative of 
sources with decreasing radial gradients in magnetic field B and electron energy distribution 
normalization factor N . Using the expressions of Marscher (1977), who considered such 
cases, the actual synchrotron turnover frequency v m is related to the intersection value v n 
by a constant factor, which in most cases is close to unity. The quantity F m is the turnover 
flux density determined by extrapolating the optically thin spectrum down to frequency 
v m . The values of u n , u m , F m , a mm , and a t hick for each source are listed in Table 8. 

For sources observed with VLBI, the angular sizes are derived from those obtained from 
the models that best fit the visibilities as described in §2.2. Though elliptical Gaussians 
are used to represent the brightness distribution of each source component, for the SSC 
calculations the visibilities can equivalently be represented by a uniform sphere with angular 
diameter 9 US = 1.8y/9 a 9b flMarscher 1987|) . The equivalent size for a spherically symmetric, 



non- uniform source is (1.8 kg) l 9 us , where kg is a parameter in the range 1-2 ( [Marscher 



1977 ). The value 9 US is used in the subsequent modeling calculations. This angular size 
value is only strictly valid if the core is known to be optically thin. In most cases, our 
spectra reveal the core to be optically thick; however, we use these optically thin values for 
consistency with the following Compton calculations (which assume the source is optically 
thin). 



3.2. Synchrotron self-Compton Calculations 

Once F m and v m are determined, these values as well as millimeter spectral indices 
and angular sizes are used with the following equation (multiplied by a constant factor; see 
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Marscher 1987] ): 



F ^ OC ^ _2 ( 2amm + 3 )_p 2(o mm +2) 



rn 



keV 



ln(v 2 /v m ) 



( 



1 + z 

5 




(3.2.1) 



Here, zaj is the upper cutoff frequency to the synchrotron spectrum and S is the Doppler 
factor corresponding to bulk relativistic motion. We use eq. (3.2.1) to predict the X-ray flux 
density at 1 keV under the initial assumption that there is no relativistic beaming (5 = 1). 
The angular sizes used in this equation should be measured at the turnover frequency. In 
general this is not the case, since VLBI observations were only undertaken at 8.4 and 22 
GHz, below v m for most sources. In the case of a turnover frequency v m > 22 GHz, the 
measured angular size is expected to follow the inequality 9{v) > 9{vm)-, and the Compton 
flux is underestimated. If the source is transparent at 22 GHz, then the measured angular 
size is acceptable. For sources with IR data, the cutoff frequency v 2 can be estimated, and 
is typically ~ 10 14 Hz. Otherwise, this frequency is assumed to be 10 14 Hz. 

In a number of cases the X-ray flux densities predicted in this way are much greater 
than those measured. If relativistic beaming is taken into account, a minimum Doppler 
factor can be derived from the measured and predicted X-ray flux densities: 



The best estimate of the Doppler factor is then used with the initial estimate of the X-ray 
flux density to generate model first and second order SSC spectra. Model spectra for 
selected FRSQ's with minimum Doppler factors obtained from the SSC calculations are 
shown in Figures 4 (0133+476) and 5 (0836+710). The parameters used to generate the 
models are listed in Table 9. The SSC calculations show that seven of sixteen sources 
observed with both VLBI and ROSAT have predicted X-ray flux densities at 1 keV greater 
than the observed values. The resolution of this discrepancy requires relativistic beaming 
of the radiation. 




(3.2.2) 
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To evaluate more clearly what these predictions mean, we can roughly estimate the 
uncertainty in the SSC predicted flux and the Doppler factor using standard propagation 
of errors techniques (i.e., |Bevington 1969| ). Though this strictly only holds for parameters 
which have normal distributions and symmetric errors, a meaningful result for the 
uncertainty can be evaluated if we treat the terms of eqs. 3.2.1 and 3.2.2 logarithmically: 

a tot = a2(J lg6 + tfrfogF + ^tfogu + ( 3 - 2 - 3 ) 

Here, a tot refers to the total logarithmic uncertainty in the predicted Compton flux, a, b, c 
refer to the powers involving spectral index which appear in eq. (3.2.1). The corresponding 
subscripts refer to angular size, maximum synchrotron flux density and the spectral 
turnover frequency. There are then other smaller or unquantifiable terms which we have 
not written down (but we discuss them below). Following this technique, and assuming 
approximately 50 % errors in the angular size, and 40% errors in the turnover flux and 
frequencies, and a spectral index of approximately 0.7, the uncertainty in log F UXtPre d is 
about 2. Reallistically, since the spectral index itself can be uncertain by a few percent, 
and since there are very likely to be correlated error terms which we haven't included, 
this number could be closer to 3. Though correlated error terms could also reduce the 
overall uncertainty, in the absence of more information we feel that it is wise to adopt the 
larger number. This number corresponds to a "worst case scenario" . As discussed earlier 
in §2.2, some of the VLBI data have angular sizes with uncertainty as small as 13%. The 
resulting uncertainty would be smaller by more than an order of magnitude in those cases 
(~ 50). However, in general, over-estimations of the Compton flux, are only significant 
if this is an over-evaluation by more than 3 orders of magnitude. By this standard, all 
of the over-estimates in Table 9, except for 0133+476, are marginal, however many of 
the underestimates are not. Using similar analysis on 5 m i n , and realizing the F vx ^ Tei j_ is 
the main contribution to the uncertainty, we calculate that the uncertainty in the log of 
the minimum Doppler factor is about 0.3. Thus, we conclude that the uncertainites are 
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too large, generally, to determine whether Compton scattering and Doppler beaming are 
necessary to reconcile predictions with the X-ray data. 

In some cases the X-ray spectral slopes predicted by the SSC mechanism disagree 
with the slopes determined from fits to the ROSAT data. This could be due to non-SSC 
contributions (e.g., from a "soft X-ray excess" component) or simply a poorly derived 
X-ray spectrum arising from the limited energy range to which the ROSAT PSPC is 
sensitive. Nonuniform models (e.g., shock waves in which the electrons suffer radiative 
losses; [Marscher fc Gear 1985|) can produce self-Compton X-ray spectra that are flatter 
than the optically thin synchrotron spectra. 

For 0836+710, the 7-ray spectrum of 1992 March is roughly in agreement with the 
predicted composite first and second order multiwaveband spectrum (Figure 5). These 
results are similar, in general, but not in detail, with the modeling results for 0836+710 
presented by Comastri et al. (1997). The model predicts that the source should have 
become faint at 7-ray energies by 1992 Oct-Nov, and indeed no detection was reported when 
the source was next observed in early 1993 ([Thompson et al. 1995|) . This illustrates the 



value of contemporaneous multiwaveband observations in testing models: source variability 
invalidates tests based on non-contemporaneous measurements. 

In addition to modeling the spectra, we use observed synchrotron parameters, such 
as F m , u m , a, z/2, and 8, in Table 10 to determine physical parameters such as iVo and B. 
These results show that the magnetic field B ranges from 1 x 10~ 2 to 140 (£ m i„/<5) G and 
iV from 9 x 10~ 13 to 2 x lO' 3 (5 / 5 min ) 2(a+2) (cgs units). 

For sources previously observed to exhibit superluminal motion, the range of possible 
values for the Doppler factor derived from the apparent transverse velocities f app can be 
compared to those derived using eq. (4.2.2). Known superluminal sources with measured 
X-ray fluxes lower than the predicted values tend to have minimum required values of S 
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that are compatable with the observed values of v app /c, as should normally be the case: 
0836+710 (v app /c = 15/1- 1 IIKrichbaum et al. 1990|| , S > 3), 4C 39.25 (v app /c = 3.5/r 1 , 
5 > 5.2 ||Marscher et al. 1991|| , and 3C 345 {v app /c = 8h~ l , 5 > 2.1; see |Unwin et aQ99j) . 
(Here we adopt the values qo=0 and Hq =100/i km s -1 Mpc -1 .) However, this need not be 
the case for all sources, since the value of 5 can be high while v app /c is low for jets pointing 
essentially along the line of sight. In such cases, the jet should not appear as extended as 
in the high v app /c sources. 



4. Statistical Analysis 

We use correlation and linear regression to test for and to model approximately the 
bivariate relationships among the data. We employ these techniques with emphasis on 
those parameters most likely to be related to X-ray emission and 7-ray emission. Since 
the SSC process predicts a specific relationship between high energy flux densities, radio 
flux density, spectral turnover frequency, and angular size, these parameters are studied 
most closely. In addition, other parameters such as the X-ray and millimeter-wave spectral 
indices, and VLBI compactness, are examined. 

All of the luminosities in this section are calculated assuming qo=0 and H —100h 
km s _1 Mpc -1 , with h arbitrarily chosen to be 0.75. 



4.1. Survival Analysis— Correlation and Regression 

Due to the presence of upper and lower limits to some of the measured or inferred 
quantities (data with limits are often referred to as "censored" data), survival analysis 
techniques are used to evaluate more accurately correlation coefficients (with their associated 
probabilities) and regression slopes fllsobe et al. 1986| ). In general, these techniques take 
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into account censored data by assuming a form for the distribution of the censored values 
relative to the measured ones. The ASURV Rev 1.1 statistical software package is used to 
calculate correlation properties for censored data Qlsobe fc Feigelson 1990|) . As pointed out 



by |Isobe et al. 1986] , spurious correlations (or spurious null correlations) can occur if these 



limits are disregarded. 

There are many correlation tests and linear regression methods. In general, since they 
are both non-parametric, the Spearman rank correlation coefficient and Kendall's r are 
best suited to test for correlation among two variables for which the nature of the parent 
distributions is unknown. The Spearman rank coefficients are provided in the summary of 
our results in Table 11. 

The Buckley- James method (BJ) is used for determining regression coefficients. This 
method is preferred over other methods, since it does not assume a Gaussian distribution 
of the residuals about the regression line. In cases for which there are limits present in 
the independent variable (e.g., the possible dependence of X-ray flux density on angular 
size), only Schmitt's regression method (SB) can be used. However, this method is not 
considered preferable for general use, since it involves arbitrary binning of the data. (For 
a more detailed discussion of the relative merits of these methods see [Isobe et al. 19861) . 



When no limits are present, the regression results are nearly identical to the results from 
standard methods (e.g., |Bevington 1969|) . 



Since several previous studies have shown a strong correlation between X-ray and radio 
luminosities for radio loud quasars (e.g., [Worrall et al. 19"87|) , we examine this relationship 



first. We consider separately the correlation between X-ray luminosity and radio core 
luminosity, X-ray and millimeter luminosity at 270 GHz, and X-ray and radio luminosity 
at spectral turnover. The analagous relationships between the flux densities as well as 
between the flux density at each waveband and redshift are studied, to ascertain whether 
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the luminosity correlations are intrinsic or spurious. 

As can be seen in Figures 6-11 the X-ray and millimeter flux densities are correlated, 
the X-ray and radio core flux densities are uncorrelated, and the X-ray and turnover flux 
densities are weakly correlated. All of the luminosity-luminosity plots in these figures show 
moderate or strong correlations. The X-ray /millimeter relationship may be stronger than 
that of the X-ray/core because the core is optically thick. Optically thick synchrotron 
emission only makes a very small contribution to the total soft radiation field which is 
upscattered to X-ray energies ( [Band fc Grindlay 1986| ). A summary of the correlation and 



regression analysis is given in Table 11. Regression results are given only if the probability 
that a correlation is not present is less than 20%. 

If the correlation between the luminosities were due solely to the common distance 
factors used to derive them, then the regression slope would be close to unity. However, the 
regression slope is less than one in all cases. Isobe et al. (1986) show through Monte Carlo 
simulations of data from flux limited samples that if the underlying physical relationship 
between the luminosities is of the form L 2 = KL 1 ^, then luminosity diagrams will reflect 
this through the regression slope, as long as upper limits are included in the analysis. The 
resulting flux-flux plots show a correlation with considerably steeper slopes. Although 
the underlying physical relationship between the luminosities in this sample is likely to 
be more complex than the one assumed by Isobe et al. (1986) their results do indicate 
that, for some forms of luminosity functions, the data will show the intrinsic correlation 
in the luminosity-luminosity plots even if the sample has a distance-luminosity bias. To 
evaluate the possibility of spurious luminosity correlations due to common redshift factors, 
we have used the partial correlation methods of Padovani (1992). This analysis calculates 
the correlation between two variables, while effectively removing the effects of other known 
variables. We find that there is still a moderate correlation between L x and both L mm and 
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L c {p s — 0.28 and 0.39 respectively) , after the correlations of each luminosity with redshift 
have been negated. However, the partial correlation between L x and L m is weak, indicating 
that the apparent relationship between these two luminosities found earlier is probably 
spurious. The results are not strongly affected if the extremely low redshift objects 3C 111 
and 3C 120 are removed from the sample (as is seen in Table 11). 

Figures 12-17 give correlation plots and histograms of various other interesting 
observational parameters. Owen et al. (1981) used the histogram of the effective millimeter 
to X-ray spectral index to show that the X-ray fluxes are closely linked to the millimeter 
fluxes. Bloom & Marscher (1991) showed that this also holds in general for the radio 
cores of quasars, even though the observations were often separated by years. Figure 15 
shows that this is also true for the present sample, with millimeter and X-ray observations 
separated by 6 months at most. Inspection of Fig. 12 and Tables 11 & 12 show that there 
is a weak anti-correlation between X-ray spectral index and both millimeter and infrared 
spectral indices. This would tend to show that a harder component to the X-ray emission 
is visible when the synchrtron spectrum is steep, and not contributing as much to the 
X-ray flux via the SSC mechanism. This finding is somewhat surprising, since we would 
expect that sources with less significant SSC emission would then be dominated by the 
softer (steeper) spectrum of an accretion disk. Though, since several of the sources with 
the steepest millimeter spectral indices are at high redshift, it could be that the softest 
X-rays from an accretion disk are redshifted out of the ROSAT band. The anti-correlations 
of the spectral indices are especially interesting considering that the fluxes in the millimeter 
and X-ray wavebands (and the luminosities) are correlated. It would be hard to argue in 
favor of the SSC mechanism (at least in a simple form) if these correlations are observed 
in the future for simultaneous data from a large set of objects. However, without detailed 
X-ray spectral analysis of individual sources (which will have to await AXAF) it is difficult 
to analyze this result much further. There is also a strong correlation between a mmx and 
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VLBI core fraction at 8.4 GHz, as seen in Figure 16. Here, core fraction is defined to be 
the fraction of total milliarcsecond-scale flux density (as determined by model fitting) that 
is in the core component. Note that there is one outlying point, which corresponds to 3C 
120. Even with 3C 120 excluded the correlation is highly significant and the regression 
relation remains similar to the earlier case. The correlation weakens at 22 GHz (Fig. 17). 
Although this may be due to the inferior calibration at high frequencies, it is also true 
that the jet (as opposed to core) components tend to be weaker, and the core stronger, at 
22 GHz than at lower frequencies. The correlation suggests that more X-rays are created 
for a given mm- wave luminosity if jet components are more prominent. That sources with 
more prominent jet components should be brighter at high energies was also suggested by 
Wehrle & Mattox (1994) in their correlation between VLBI structure and 7-ray detections, 
although their result was of marginal statistical significance. It is certainly not expected if 
the bulk of the high-energy emission comes from the quiescent jet, as in the Ghisselini et 
al. (1985) model. 

There are various explanations for this correlation. If the bulk Lorentz factors for 
each source are similar, then the sources with smaller angles between the jet axis and the 
line of sight will have larger Doppler beaming factors. Such sources will be more severely 
foreshortened as viewed by the observer, and thus the jet "knot" components will not be as 
easily detected (because of blending with the core) as for jets that are less highly beamed. 
This effect will also tend to decrease the SSC X-ray flux for a given millimeter flux (all 
other quantities, such as turnover frequency, being roughly equal; see eq. [3.2.1]). Recent 
VLBI and X-ray observations of 3C 345 ( [Unwin et al. 1997| ) indicate that there are states 
in which the X-ray emission emanates from the jet and not the nucleus, as one might 
normally expect. At these times the relative contribution of radio flux from the extended 
jet is much higher. However, we note that a mmx was uncorrelated with the core fraction as 
this particular flare progressed in 3C 345. 
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Another possible explanation is that the sources with less prominent jets have wider 
opening angles, so knots expand rapidly, which would decrease their SSC emission. If this 
is the case, such sources should be found to be more highly variable (at least in the decay 
stages of flares) than those with prominent jets. It may also be that the sources with 
strong jets are in an enhanced state (e.g., from shock compression), which would increase 
the ratio of X-ray to millimeter flux as compared to the other sources. However, for several 
individual sources (NRAO 140, 0836+710, and 4C 39.25), this ratio seems to be preserved 
during flares. As mentioned in §2, long term X-ray variability appears to be related to 8.4 
GHz VLBI structure as well, though this information is probably too limited. In the Melia 
& Konigl (1989) model, the X-rays arise from inverse Compton scattering of seed photons 
from an accretion disk by relativistic electrons deep in the jet. Because of the geometry of 
the scattering, the X-ray emission is lower along the jet axis than at a small angle (roughly 
the inverse of the Lorentz factor of the jet); viewed from this angle, the jet components 
would have maximum visibility. 

There are no significant correlations between X-ray flux density and the other 
parameters used to predict the X-ray flux. This is somewhat surprising since one might 
expect strong correlations from relations such as eq. (3.2.1). A lack of correlation between 
X-ray flux density and angular size can, however, be understood since many of the angular 
sizes are actually upper limits, and other angular sizes may represent the blending of two 
or more components. Futhermore, the actual angular sizes at the turnover frequency are 
likely to be smaller than those measured at 22 GHz (see §3). The weak correlation between 
X-ray flux density and turnover frequency can be understood since many of the turnover 
frequencies are determined with limited data and are approximate. This uncertainty would 
then also be exacerbated by an inaccurately determined turnover flux density. Of course, 
the simplest explanation for all of these non-correlations is that the X-ray flux is not SSC 
in origin, and thus none of the correlations would be strong. If this is the case, the actual 
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Doppler factors must exceed the minimum values obtained from the SSC analysis (Table 9). 



4.2. Distributions 

The histograms of the observed properties of FSRQ's show which values are most 
common among the objects. These reveal both selection effects and the physical aspects of 
the objects, as is demonstrated below. The distribution of turnover frequencies (Fig. 13) 
shows a very broad peak at roughly 50 GHz, similar to that found by Gear et al. (1994). 
The distribution cuts off well below 150 GHz at high frequencies. If the flat spectra are 
caused by a jet with radial gradients in magnetic field and electron density ([Marscher 1977 



discusses the spherical case; see [Marscher 1980) ; [Konigl 1981| for application to jets), then 



the angular size of the emission region, observed in the partially opaque portion of the 
spectrum, decreases inversely (the precise dependence is a function of the gradients in B 
and N) with frequency up to the spectral turnover. Thus, the spectral turnover frequencies 
are indicators of the innermost radius of the source. That the spectral turnovers occur 
at frequencies ~ 50 GHz, rather than up to an order of magnitude higher (as has been 
observed for some BL Lac objects and flaring quasars; see [Brown et al. 1989] ; |Gear et aT| 



1994|) indicates that the smallest radius is farther out from the central engine than is implied 
for these other sources. The relationship 6 m oc 622(v/2'2GHz)~ k2 is used to approximate 
the size of the inner core for the sources in our sample, where #22 is the angular size 
measured at 22 GHz. If that size is an upper limit, a value equal to half of the upper limit 
is used. The exponent k 2 is a constant of order unity that depends on the optically thick 
spectral index and is enumerated in Marscher et al. (1977). A histogram of this estimated 
"innermost radius" is provided in Figure 14. The peak in the distribution of turnover 
frequencies (see Fig. 13) is an indicator of the peak in the histogram of the inner radius 
at ~ 0.6 parsecs. Though, considering our previous discussions (§2.2; §3.1-3.2) regarding 
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the uncertainty of the angular size, assuming an optically thin/thick source, and possible 
blending of components, we warn that this number should be seen as a rough estimate of 
where a physical change in the jet should take place and not as a canonical number. The 
physical meaning of this inner radius is discussed in §5 below. 

The distribution of millimeter spectral indices (see Table 12) has a sharp peak at ~ 
0.7, and is clearly skewed toward steeper spectral indices. 



5. Discussion 

We have shown that the high energy properties of FSRQ's are related to their radio to 
infrared properties. In particular, the mm-X-ray luminosity correlation is strong, as found 
in previous studies. In addition, we find an anti-correlation between X-ray flux and VLBI 
core dominance, which suggests that the knots which comprise the observed radio jets in 
FSRQ's contribute substantially to the X-ray emission. 

That the distribution of spectral indices peaks near a mm ~ 0.7 indicates that most of 
the sources contain a distribution of relativistic electrons with a canonical E~ 2A spectrum. 
However, several sources have much steeper millimeter-wave spectral indices, as expected 
for radiative energy losses, or emission from an inhomogenous jet ( Marscher 1980f ). This 



point is further borne out by the infrared observations. Of the fifteen sources for which we 
have obtained infrared data (reported here in Table 6 and in |Bloom et al. 199~3| ), six show 



infrared spectral indices that are 0.5 to 1.0 steeper than the mm-wave spectral indices, 
which is again consistent with the basic predictions for radiative energy losses. However, 
seven sources have infrared spectral indices that are flatter by similar amounts, probably 
because of another (thermal?) emission component becoming prominent toward higher 
frequencies. 
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The combined peaks in the distributions of spectral turnover frequency and angular 
size are shown above to indicate a characterictic size of about 0.6 pc for the radio cores of 
the FSRQ's observed with VLBI. In an inhomogeneous jet model with radially decreasing 
magnetic field and electron energy normalization, the synchrotron turnover frequency is 
predicted to increase with decreasing radius (Marscher 1980| ; |Ghisselini et al. 1985|) . Hence, 
if jets extended down to some arbitrarily small radius, this would result in much higher 
observed turnover frequencies. Thus, there is clearly some physical inner cutoff to the radio 
jet. The exact nature of this is unknown. There could be a change in geometry at this 
radius. That is, the jet could make a transition from a conical geometry to one which 
has a weaker dependence on radius, such as a paraboloid (e.g. |Marscher 1980| ). Magnetic 



field and electron density, and therefore the turnover frequency, would then depend less 
strongly on radius. This transition in geometry could be caused by a change in the external 
pressure on the jet or the cooling of the protons to nonrelativistic mean internal energies. 
A change in Doppler factor along the jet, either from a change in jet speed or orientation, 
could also have the effect of creating a cutoff in the emission since the observed flux from 
the jet follows F u oc 5 3+otr . Millimeter- wave VLBI observations are important to study this 
transition region. 

The lower redshift objects appear to have distinct X-ray properties. For example, 
both 3C 111 and 3C 120 have very steep X-ray spectra (a x > 2) and the X-rays could 
have a thermal origin ( Permer et al. 1993| ; |Grandi et al. 1997|) . This steep soft X-ray 
spectral component would lie at frequencies below the ROSAT PSPC response at the 
higher redshifts of the quasars in the sample. Indeed, at higher X-ray energies both radio 
galaxies have flatter spectra, similar to quasars (Petre et al. 1984] ) . Leach et al. (1995) 
have shown that the spectral slope of the low- redshift quasar 3C 273 is —2.7 in the low 
energy band (0.1-0.3 keV) of the ROSAT PSPC and -0.5 in the hard band (1.5-2.4 keV). 
Also, as discussed earlier, 0736+017, at a redshift of 0.19, has an X-ray spectral index ^2 if 
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photoelectric absorption by Galactic molecular gas is taken into account. 

Three sources, 0836+710, NRAO 140, and 4C 39.25, clearly show radio to millimeter 
variability that is directly proportional the X-ray variability. As discussed by Marscher 
(1988), this would occur if the quantity N x is nearly constant with time, with N varying 
inversely as the line of sight distance through the source x. However, an increase in the 
Doppler factor would also have a similar signature (see [Marscher et al. 1992| and below). It 



is further likely, from the general considerations of Bloom & Marscher (1996) and the SSC 
calculations in §3, that the high states of radio to X-ray (and in some cases 7-ray) emission 
are associated with higher high-energy cutoffs to the electron distribution. This increases 
the X-ray flux density by another factor of 21n7 maa; . This would predict that if the 7-rays 
are from first order SSC, then the X-ray flux density should increase by a slightly larger 
factor than the synchrotron flux density. In cases for which the entire spectrum increases 
by a fixed amount, the increase can be attributed to a change in the Doppler factor. 

Models that involve the scattering of external radiation from the relativistic electrons 
in jets (e.g., Melia and Konigl 1989; Dermer & Schlickeiser 1993) would also explain this 
linear relation between millimeter and X-ray flux if changes in Nq were solely responsible 
for the variations in the synchrotron and scattered fluxes. This would be the case if density 



fluctuations occurred while the magnetic field remained constant (see Bloom fc Marscher 
19961) . 



Of the six sources in our sample that were detected in high energy 7-rays, two were in 
a higher state at mm or cm wavelengths at the time of the detection. 0836+710 showed 
a factor of two increase in the millimeter and X-ray flux densities when it was detected 
in 7-rays relative to earlier times. Though 0804+499 was only observed at millimeter 
wavelengths several months before and after the 7-ray detection, cm-wave flux densities 
clearly increased by a factor of 3 between August 1991 and December 1992. The other 
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sources (0917+449, 0954+556, 1611+343, 1633+382) were not well observed during the 
time of the EGRET observations. 

The physical relationships among the various parameters are often complicated and 
hidden by selection effects. We have used statistical tests to attempt to reveal underlying 
relationships. We find that the X-ray luminosity is correlated with the radio core and 
millimeter-wave luminosities. However, the strong correlation is probably due in part to 
luminosity-redshift bias. 

There are many aspects of this work that should be continued to gain a better 
understanding of the relationship between the radio to IR and X-ray to 7-ray properties of 
FRSQ's. Because many sources in this sample exhibit spectral turnovers near 40 GHz, and 
these turnovers indicate that the emission originates from the innermost region of the jet, 
they should be observed with VLBI at 43 GHz (VLBA) and in the 1-3 millimeter range 
(various arrays) to try to resolve these inner regions and gain a better understanding of the 
physical conditions at the base of the radio jet. In addition, the angular sizes derived from 
these VLBI images can be used for more accurate Compton calculations in conjunction 
with X-ray and 7-ray observations. Simultaneous multiwaveband and VLBI observations 
would relate the variations at different frequencies both to each other and to changes in 
the structure of the jet. Time-dependent models of SSC emission from relativistic jets 
should be developed to replace the crude uniform spherical source approximation used in 
this study. Finally, complete samples chosen using various well defined criteria should be 
studied thoroughly to establish statistical relationships as well as to elucidate the effects 
of those selection criteria. Together, future studies such as these should lead to a more 
complete understanding of the nature of the nonthermal emission in FRSQ's. 
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FIGURE CAPTIONS 

Fig. I. — VLBI images at 8.4 and 22 GHz from the Global VLBI Network. Parameters of 
each map are listed in Table 3. The restoring beam is shown in the lower left corner of each 
separate plot. 

Fig. 2.— VLBA images at 22 GHz 

Fig. 3. — Determination of the turnover frequency from VLBI and total flux density data. 
The dotted lines show the frequency v n at which the extrapolation of the partially opaque and 
optically thick spectra meet. The solid curve is the synchrotron spectrum for a homogeneous 
sphere with spectral turnover at frequency v m as discussed in the text. 

Fig. 4. — The synchrotron spectrum (fit by a self-absorbed, uniform source model) and 
predicted first order (middle curve) and second order (rightmost curve) spectra for 0133+476, 
plotted against the contemporaneous data. 

Fig. 5. — Similar to Figure 4, except for the quasar 0836+710. The EGRET gamma-ray 
data from Thompson et al. (1993) is plotted to the right on the 1992 March spectrum. 
Optical data are from von Linde et al. (1993). The large uncertainty in the optical point 
corresponds to the range of fluxes during a short term flare. 

Fig. 6. — Correlation plot between VLBI core flux density at 22 GHz and X-ray flux density 
at 1 keV. 

Fig. 7. — VLBI core luminosity at 22 GHz vs. X-ray luminosity at 1 keV. The solid line 
corresponds to a linear regression. 

Fig. 8. — Millimeter-wave flux density (273 GHz) vs. X-ray flux density at 1 keV. 
Fig. 9. — Millimeter-wave luminosity (273 GHz) vs. X-ray luminosity at 1 keV. 
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Fig. 10. — Flux density at the turnover frequency v m vs. X-ray flux density at 1 keV. 

Fig. 11. — Luminosity at the turnover frequency u m vs. X-ray luminosity at 1 keV. The 
solid line is a linear regression. 

Fig. 12. — X-ray spectral index vs. millimeter- wave spectral index. 
Fig. 13. — Histogram of turnover frequencies v m . 

Fig. 14. — Histogram of innermost radius measured by our 22GHz VLBI observations. 

Fig. 15. — Histogram of X-ray to millimeter-wave spectral index, defined as a mmx = 
\og(F x /F mm )/\og(2A x 10 17 /2.7 x 10 11 ) 

Fig. 16. — X-ray to millimeter-wave spectral index a mmx vs. radio core fraction at 8.4 GHz. 
Fig. 17. — X-ray to millimeter- wave spectral index a mmx vs. radio core fraction at 22 GHz. 
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TABLE 11 

Correlation and Regression Parameters 
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